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1.0  INTRODUCTION 


Til e  objective  of  this  research  effort  was  to  determine 
the  feasibility  of  making  significant  improvements  in  seismic 
depth  phase  detection  using  statistical  signal  processing  tech¬ 
niques  developed  at  PNSCO,  INC.  Applications  of  these  tech¬ 
niques  together  with  additional  techniques  developed  in  the 
course  of  the  investigation  indicate  significant  depth  phase 
detection  improvements  are  possible. 

The  primary  interest  in  detecting  seismic  depth  phases  is 
that  good  estimates  of  event  depths  can  be  obtained  from  the 
(pP-P)  and  (sP-P) arrival  time  differences.  This  type  of  depth 
determination  is  preferred  over  all  other  depth  estimation  pro¬ 
cedures  . 

At  present,  depth  phase  detection  is  usually  performed  by 
visual  identification  from  the  seismogram  and  is  not  successful 
for  many  events.  This  type  of  detection  uses  only  a  very  limited 
portion  of  the  depth  phase  information  contained  in  the  seismo¬ 
gram  and  is  successful  only  when  depth  phases  arc  particularly 
strong.  liven  when  depth  phases  are  visually  detected  there  can 
be  substantial  errors  in  the  delay  time  estimates  since  they  are 
based  on  the  identification  of  P,  pP,  and  sP  first  movements. 

A  procedure  which  utilizes  a  larger  portion  of  the  P,  pP  and  sP 
waveforms  would  provide  a  more  precise  estimate. 
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To  date,  depth  phase  detection  based  on  techniques  such 
as  the  auto- corre 1  at  ion ,  cepstrum,  etc.  have  also  had  limited 
success.  Although  these  techniques  use  more  of  the  informa¬ 
tion  contained  in  the  seismogram,  detection  commonly  deterior¬ 
ates  rather  than  improves  when  more  of  the  coda  is  used  in  the 
analysis. 

The  statistical  processing  techniques  which  we  have  de¬ 
veloped  are  based  on  simple  physical  models  and  insure  construc¬ 
tive  use  of  the  dep th  phase  information  contained  in  the  entire 
seismogram.  These  new  techniques  are: 

1.  Stochastic  stacking  -  this  technique  allows 
for  deptli  phase  delay  time  variations  along 
the  coda  which  arc  associated  with  the  later 
seismic  arrivals. 

2.  Phase  stacking  -  the  consistency  of  the  phase 
difference  between  the  direct  and  reflected 
waves  can  be  used  to  enhance  depth  phase  de¬ 
tection. 

Recordings  of  a  seismic  event  were  analyzed  at  six  sta¬ 
tions  using  the  new  statistical  techniques  as  well  as  several 
of  the  conventional  techniques.  for  the  event  analyzed,  the 
depth  phases  were  not  visually  apparent.  The  conventional  sig¬ 
nal  processing  techniques  detected  the  pi1  at  one  of  the  six 
stations  whereas  the  new  techniques  detected  both  the  pi’  and  sP 
at  one  station,  and  the  pP  at  two  others. 
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It  was  also  determined  that  by  stochastically  averaging 
the  data  from  several  stations  depth  phases  could  be  detected 
that  were  not  apparent  from  analysis  of  the  data  from  any  of 
the  individual  stations. 

The  results  based  on  this  work  indicate  that  these  new 
techniques  could  greatly  increase  the  frequency  and  accuracy 
of  the  detection  of  seismic  depth  phase  delay  times. 


2.0  Sl'lSMJC  DliPTll  PI  I  ASH  DLTFCTION 


2  •  1  Simplified  Geophysical  Mode ] 

The  simplified  geophysical  model  illustrated  in  Figure  1, 
can  help  indicate  some  of  the  problems  associated  with  using 
tlie  entire  record  for  seismic  depth  phase  detection. 

In  general,  a  seismogram  is  constructed  of  seismic  arri- 
\uls  taking  many  different  paths  from  source  to  receiver.  The 
first  of  these  seismic  arrivals  is  the  direct  P  wave  followed 
by  its  associated  reflected  phases  pP  and  sP,  which  arrive  a- 
long  minimum  time  paths  very  similar  to  that  of  the  P  wave  for 
event  depths  of  "50  km  or  less.  It  is  the  arrival  time  differ¬ 
ences  of  these  first  three  arrivals  which  arc  most  useful  for 
depth  determination.  (In  Figure  2  arc  examples  of  seismic  re¬ 
cordings  for  which  the  P,  pP  and  sP  phases  arc  easily  identified.) 
At  some  later  point  in  time  seismic  waves  which  reflect  off  the 
core  arrive  as  Pel’,  pPcP  and  sPcP  phases.  Later  yet  (depending 
on  the  event  location)  PP,  pPP  and  sPP  phases  will  arrive.  The 
point  here  is  that  the  seismogram  is  constructed  of  triplets  of 
arrivals  from  each  of  the  many  seismic  paths  from  source  to  re¬ 
ceiver.  In  general  the  time  delays  within  these  triplets  will 
differ  since  the  relative  geometry  of  the  paths  within  each  tri¬ 
plet  is  different.  The  consequence  of  this  is  that  the  depth 
phase  delay  time  will  vary  as  a  function  of  position  along  the 
coda.  As  an  example,  for  a  -10  km  event  at  45°,  the  pP-P  delay 
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is  'll  sec.  whereas  the  pPP-PP  delay  is  -10  sec.  Adjusting 
for  these  variations  by  use  of  travel  time  tables  is  not  an 
obvious  procedure  since  the  relative  strength  of  the  various 
overlapping  phases  and  the  event  depth  must  be  approximately 
known . 


To  sec  how  such  delay  time  variations  can  deteriorate 
depth  phase  detection,  consider  the  cepstrum  analysis  of  a 
synthetic  seismogram.  The  cepstrum  is  a  near  optimum  echo 
detector  for  Gaussian  echoed  signals  in  Gaussian  noise. 

In  figure  3  arc  the  cepstrums  of  synthetic  seismograms  for 
events  at  40  and  15  km  depth;  the  horizontal  axis  is  time  in 
seconds.  The  synthetic  seismograms  were  generated  by  passing 
white  noise  through  a  recursive  digital  filter  having  a  band 
pass  typical  for  seismic  arrivals  and  adding  this  signal  to 
itself  at  delays  and  (corresponding  to  pP-P  and  sP-P) . 

The  recursive  digital  filter  was  designed  from  the  following 
z- transform  of  a  resonator  with  poles  at  z  =  re-  >rT  and  a  zero 
at  q; 


H  ( z)  = 


1  -  q  z 


-1 


(1  -  2r  (cosuj^.T)  z" 1  +  r2z"2) 


Foi  q-1  this  gives  zero  gain  at  w-0  and  a  resonant  re¬ 
sponse  with  ur  determining  the  resonant  frequency  and  r 
relating  to  the  Q  of  the  response.  We  used  q=l,  r=.9  and 
wr=2"  for  the  generation  of  the  synthetic  data  with  T  being 
the  inverse  of  the  sample  rate. 


In  this  simulation  the  delay  times  and  t ^  were 
constant  throughout  the  seismogram  and  clear  ccpstral  peaks 
indicate  the  delays  ^  and  t 2  for  the  40  km  ease.  Peaks 
Vv i  1 1  also  appear  at  many  of  the  possible  combinations  of 
arrival  times  (nt^  +  m t 2 )  where  n  and  m  are  integers.  In 
the  15  km  example  it  is  apparent  that  interpretation  of 
these  results  should  be  done  using  the  entire  ccpstrum 
pattern . 


If  t ^  and  t 2  vary  along  the  seismogram  by  times  comparable 
to  the  Jin  If  widths  of  the  ccpstral  peaks  or  greater,  the  pro¬ 
minence  of  these  peaks  would  be  substantially  reduced  since  their 
amplitudes  would  be  distributed  over  a  range  of  delay  time  values. 
To  compensate  for  this  type  of  reduction  in  ccpstral  peak  detec¬ 
tion  a  stochastic  stacking  technique  is  used. 


2 . 2  Stochastic  Stacking 


A  stochastic  stacking  technique  can  be  used  to  increase 
the  detectability  of  a  peak  whose  position  can  vary  unprc- 
dictably  within  a  certain  limited  range.  To  see  this,  first 
consider  a  set  of  cepstrums  calculated  from  consecutive  time 
segments  of  a  seismogram.  If  the  depth  phase  delay  time  were 
the  same  in  each  segment,  then  by  stacking  (adding)  these  cep¬ 
strums,  the  amplitudes  of  the  stationary  peaks  would  add  where¬ 
as  the  non-stationary  peaks  arising  from  origins  other  than  the 
depth  phase  delay  times  would  average  to  some  lesser  amplitude 
level.  This  is  also  true  for  the  same  calculation  of  a  single 
ccpstrum  of  the  entire  seismogram. 


By  redefining  the  N  ccpstrum  values  x^,  i  =  1,N  to  be 


yi  *  MAX  (x. ,  j  =  i  -  a/2,  i  =  i  +  A/2) 

foi  cacn  point  i,  and  then  by  adding  the  y  arrays  one  can  in- 
cieasc  detectability  for  cases  in  which  the  peak  of  interest 
mo\ es  unp redi c t ab 1 y  within  a  window  A.  figure  4  illustrates 
this  technique  using  synthetic  data.  bach  of  the  top  six  arrays 
arc  constructed  of  random  numbers  having  values  between  0  and  1. 
A  peak  (indicated  by  an  arrow  in  Figure  4),  defined  by  three 
points  having  values  .8,  1.2  and  .8,  is  added  to  each  array,  but 
shifted  on  additional  point  to  the  right  in  each  array.  The 
lcsult  of  adding  these  arrays  (straight  stacking)  is  shown  and 
it  is  seen  that  a  clear  detection  of  the  peak  is  not  achieved 
bv  ordinary  stacking.  However,  if  before  adding,  one  redefines 
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Straight  Stack 


Stochastic  Stack 


Demonstration  c*f  stochastic 
stacking  using  svnthrtic  data 


Figure  4 


a-V-A 


these  arrays  as  has  been  outlined,  the  results  show  that  a 
dramatic  improvement  in  peal;  detection  is  achieved  by  using 
this  stochastic  stacking  procedure.  (In  Figure  A1  of  the 
appendix,  one  can  compare  this  stochastic  stacking  procedure 
with  that  of  low  pass  filtering  the  arrays  with  a  cutoff  fre¬ 
quency  equal  to  the  inverse  of  the  stochastic  time  window  A. 

It  is  seen  that  stochastic  stacking  and  linear  filtering  have 
different  effects  on  data.) 

Thus,  the  stochastic  stacking  technique  should  bo  effcc- 
ti\c  in  detecting  depth  phase  delay  times  since  these  delay 
times  vary  as  a  function  of  position  along  the  seismogram. 

2 . 3  Phasor  Stacking 

I  up ] css i vc  i mpi ovements  in  seismic  depth  phase  d  tcction, 
m  addition  to  those  gained  by  stochastic  stacking,  have  been 
obtained  b)  using  information  concerning  the  phase  difference 
between  the  direct  and  reflected  waves.  To  sec  how  this  can 
be  achieicd  consider  the  steps  involved  in  calculating  the  ccp- 
strun. 

Consider  the  received  signal  F(t)  to  be  the  sum  of  the 
direct  wave  f(t)  and  a  single  echo  of  relative  amplitude  a, 
phase  0,  and  echo  delay  time  t.  F(t)  can  be  written 

F(t)  =  f(t)  +  a  (f(t  -  t)  cos  0  +  f{[  (t  -  t)  sin  0] 
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Here  fn  represents  the  Hilbert  transformation  of  f  which  cor¬ 
responds  to  shifting  each  Fourier  component  of  f  by  it/2.  Tbc 
bracketed  expression  then  represents  the  signal  f(t)  having 
e; eh  of  its  Fourier  components  shifted  by  0.  This  phase  shift 
0  is  the  phase  difference  between  the  direct  and  reflected  sig 
nal.  lor  example,  if  the  echo  differed  from  the  direct  signal 
by  only  a  change  in  sign,  0  would  equal  tt  . 

The  power  spectrum  of  F(t)  is, 

P(w)  -  p(m)  [1  +  +  La  cos  (cot  -  0)] 


for  w  >  0  and  with  p(u)  being  the  power  spectrum  of  f (t) .  As 
can  be  seen  from  the  expression  for  P(u)  ,  the  power  spectrum 
p(w)  is  modulated  by  a  cosinusoidal  function  having  frequency 
t  and  phase  0.  Taking  the  complex  spectrum  of  P(co)  one  would 
then  obtain  a  peak  at  lag  t  which  would  have  phase  6.  Thcrc- 
foic  one  can  assign  a  phase  G  to  each  point  of  the  ccpstrum 
which  is  calculated  by  taking  the  power  spectrum  of  log  P(io)  . 

To  see  how  this  phase  information  can  aid  in  depth  phase 
detection  consider  events  at  ~50  km  or  less.  For  such  events 
the  direct  and  reflected  phases  (.e.g.,  PeP,  pPeP,  sPeP)  travel 
similar  paths  once  outside  the  vicinity  of  the  event.  It  is 
then  renso)  able  to  assume  that  the  direct  and  reflected  phases 
would  undergo  similar  reflections,  distortions,  etc.  once  away 
from  the  source  region  and  whatever  phase  relationship  they  had 
emerging  fion  the  source  region  would  tend  to  be  preserved  in 
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route  to  the  receiver.  If  the  different  reflected  phases  had 
undci gone  similar  phase  shifts  relative  to  the  direct  waves, 
the  phase  of  the  depth  phase  ccpstral  peaks,  calculated  from 
different  portions  of  the  seismogram,  should  tend  to  be  con¬ 
stant  within  a  possible  tt  shiftf^  Spurious  ccpstral  peaks  not 
related  to  the  echo  delay  times  would  have  a  more  random  phase 
as  a  function  of  which  portion  of  the  seismogram  is  analyzed. 


These  assumptions  are  really  just  generalizations  of  ideas 
normally  assumed  about  the  direct  and  reflected  waves.  It  is 
common  to  consider  the  echo  to  have  a  waveform  similar  to  the 
diicct  wa\ c  but  having  the  same  or  different  sign  depending  on 
the  angle  of  reflection  and  on  the  signs  of  the  parts  of  the 
radiation  pattern  that  radiate  the  direct  and  the  reflected 
waves.  This  is  equivalent  to  saying  there  is  a  0  or  tt  phase 
difference  between  the  direct  and  reflected  waves.  What  wre  have 
done  hci c  is  to  generalize  this  idea  to  include  all  phase  shifts 
from  0  to  2tt  rather  than  just  0  and  tt  and  to  use  the  consistency 
of  0  ±  tt  throughout  the  seismogram  to  aid  the  depth  phase  detec¬ 
tion. 


This  phase  consistency  can  be  utilized  by  stacking  ceps t rum 
pirns oi  ai i a) s  and  computing  the  amplitudes  of  the  phasor  sums, 
by  stacking  those  phvsor  arrays  ccpstral  peaks  which  have  a  con¬ 
sistent  phase  0  ±  tt  throughout  the  stack  can  become  dominant 
even  though  their  amplitude  sums  arc  not.  This  effect  will  be 
demonstrated  in  the  analysis  of  the  data.  We  also  note  that  the 
dcgi ee  of  phase  selectivity  can  be  adjusted  by  the  manner  in 
which  the  phasors  are  summed. 
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APPLICATION  OF  TECHNIQUES  TO  A  SEISMIC  EVENT 


Analysis  was  performed  on  seismic  data  recorded  at  six 
stations  for  a  single  event  using  these  new  depth  phase  de¬ 
tection  techniques  as  well  as  several  of  the  conventional 
analysis  techniques.  The  event  chc  cn  was  the  Andreanof  Is¬ 
land  Earthquake  of  22/11/65  (origin  2025  31.1,  51.4  N, 

179.7  W,  m^  =  5.9,  sample  rate  =  20  sps)  .  The  map  of  Figure 
5  indicates  the  event  location  and  the  six  LRSM  stations  from 
which  the  seismic  recordings  were  analyzed. 

CSGS  reported  this  event  to  have  a  depth  of  40  km  from 
location  fit  estimates  and  the  results  of  our  analysis  agree 
with  this  depth  estimate.  However,  the  analysis  by  SDL 
("The  Long  Shot  Experiment") , SDL  Report  234,  Vol.  II, 
page  51)  concludes  that,  "...  The  visual  analysis  of  the 
LRSM  records  across  North  America  for  the  22  November  1965 
event  revealed  a  depth  phase  pP  which  was  consistent  from 
station  to  station  and  which  resulted  in  a  depth  calcula¬ 
tion  of  approximately  20  km,  ...".  If  one  inspects  the  re¬ 
cords  used  for  this  analysis  (Fig.  5a)  one  sees  that  this 
is  not  at  all  an  obvious  conclusion.  In  Figure  5a  the  dotted 
lines  indicate  the  approximate  positions  for  pP  for  an  event 
at  20  km  and  pP  and  sP  for  an  event  at  40  km.  From  these 
records  it  would  be  difficult  to  make  conclusions  concerning 
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any  of  the  depths.  However,  if  one  considers  the  Bolivian 
Lvent  (Figure  2)  to  be  a  good  examj le  of  the  wave  forms 
received  when  the  depth  phases  are  clearly  identifiable, 
then  the  records  of  I1V-MA  and  MN-NV  (Figure  5a)  would  be 
the  most  suggestive  indication  of  depth  phases  and  these 
would  be  in  agreement  with  the  40  km  depth  estimate. 

The  seismograms  we  analyzed  are  plotted  in  Figure  6 
and  show  a  good  signal  to  noise  ratio  eventhough  the  depth 
phases  are  not  visually  apparent. 

Analysis  was  done  using  the  following  methods: 

1.  Cepstrum  of  the  direct  arrival  portion  of  the 
seismogram  (1st  -  25.6  sec.) 

2.  Cepstrum  of  a  ~4  min.  portion  of  the  seismogram 

3.  Stacked  ccpstrums  (9  -  25.6  second  samples) 

4.  Stochastic  stacked  ccpstrums  (9  -  25.6  sec. 
samples) 

5.  Stochastic  cepstrum  phasor  stack  (9  -  25.6  second 
samples) 

Additional  methods  used  for  the  comparison  of  low  pass 
filtering  to  stochastic  windowing  were: 

6.  Stacked  auto -correlat ion  functions  of  low  pass 
filtered  data  (9  -  25.6  second  samples) 
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7.  Stacked  ccpstruui  of  low  pass  filtered  data 
(D  -  25.6  second  samples) 

The  following  steps  used  in  the  calculation  of  the  ccp- 
strum  were  optimized  using  the  synthetic  data: 

•  Select  N  sampled  (sample  rate  =  20  pts/scc) 
amplitude  values  from  the  portion  of  a 
seismogram 

XM(n) ,  n=0,  1 ,  . . . ,  N-l 
and  add  N  zeroes  to  interpolate  the  spectrum 

•  Calculate  the  amplitude  spectrum  for  positive 
frequcnc ics 


amU)  = 


1 

2.N 


2K  - 1 

I 

n=0 


XM(n)e 


2minj/2N 


3=0,1, 


N-l 


with  the  Nyquist  frequency  =  10  Hz  and  i  =  (-l)1^. 

•  Retain  only  the  lower  quarter  of  frequencies 
of  the  amplitude  spectrum  since  there  is  very 
little  energy  at  frequencies  above  2.5  Hz.  This 
leaves  the  array 


C  j ) »  j  -  o ,  1 , 


•  •  •  > 


N/4-1) 


Remove  the  mean  and  apply  a  cosinusoidal  taper 
to  the  first  10u  and  last  20%  of  the  array 
giving  the  modified  array 

* 

(Am  (j)  ,  j  =  0,  1,  .  .  .  ,  N/4-1) 

The  2  0°j  taper  on  the  higher  frequencies  was  us<  d 
to  increase  the  de-emphasis  of  the  higher  fre¬ 
quencies  and  was  based  on  indications  from  other 

works  that  the  depth  phase  spectrum  has  a  greater 
similarity  to  the  direct  phase  spectrum  below 
~  1 . 5  Hz. 


Add  N/4  zeroes  to  interpolate  the  cepstrum  giving 
the  array 


(Am  (J)  >  j  =  0,  1,  .  .  .  ,  N/2-1) 

At  this  stage  one  would  take  the  log  of  this 
A^  array  to  obtain  the  log  cepstrum;  but  our 
results  show  that  by  not  using  the  log  better 
re  ults  were  obtained  and  the  results  presented 
in  this  report  were  all  obtained  without  use  of 
the  log. 


•  Calculate  the  Fourier  Transform  of  the  A.,  array. 

M  1 
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1 

WfJ 


FM(k) 


N/2-1 

J  Am  (j)e-2«jk/(N/2) 

3-0 


where  (F^(k) ,  k  =  0,  1,  .  N/4-1)  are 

complex  numbers  representing  the  positive  fre¬ 
quency  spectrum  of  A„.  One  can  now  obtain  the 
amplitude  and  phase  of  each  cepstrum  point  k. 

The  array 

(  I Fjyj (k)  |  ,  k  =  0,  1,  ...»  N/4-1)  is  what  we 

refer  to  as  the  cepstrum  amplitude  and  the  com¬ 
plex  numbers  F^(k)  ate  referred  to  as  cepstrum 
phasors . 

•  To  calculate  the  stochastic  cepstrum  stack,  one 
then  calculates 

|FM(k)  |  =  Max  (  |FM(j)|,  j  =  k-A/2,  k+A/2) 

for  each  section  M  and  then  sums  over  the  number 
of  sections  (NS)  used  from  a  given  seismogram 
(A  is  the  stochastic  window  width)  giving 
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NS 

coo  -  Y,  iJMCk)  I  *km 

M=1 

where  C(k)  is  the  stochastic  ccpstrum  stack  and 
WM  is  a  weighting  factor,  chosen  such  that  the 
mean  amplitude  of  each  | FM (k) |  array  is  equal. 

•  To  calculate  the  stochastic  ccpstrum  phasor  stack, 
each  F^(k)  is  replaced  by  F^(n)  where  n  is  the  in¬ 
dex  of  the  Max  value  of  |P^(j)|in  the  interval 
j=k-A/2,  k+A/2.  Then  the  stochastic  ccpstrum 
phasor  stack  is  defined  by 


CP (k) = 


NS 

2 

M=1 


w, 


M 


fmoo 


3.1  Comparison  of  the  New  Depth  Phase  Detection  Techniques 
with  the  Conventional  Methods 


We  will  now  compare  the  effectiveness  of  the  various 
techniques  in  detecting  depth  phases  from  recordings  for 
which  they  arc  not  visually  apparent.  The  data  from  all 
six  stations  was  processed  in  the  identical  fashion.  In 
Figure  7  arc  the  analysis  results  for  the  station  MN-NV. 

At  the  top  of  the  figure  is  a  plot  of  the  first  portion  of 
the  seismogram.  The  next  three  plots  arc  the  ccpstrums 
calculated  in  the  three  conventional  ways  indicated.  The 
two  vertical  dotted  lines  mark  the  expected  delay  times 
for  the  pP  and  sP  phases  for  a  39  km  deep  event ,  the  depth 
we  determined  this  event  to  be.  For  an  event  at  -60°  these 
delay  times  would  be  ~1  1  . 1  and  "'15.8  respect  ivcly . 
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In  interpreting  these  ceps tj  inns  one  seeks  dominant  peaks 
at  the  appropriate  delay  times.  In  the  ccpstrum  of  the  1st  - 
25.6  seconds  of  data,  a  peak  appears  at  the  sP  time  delay;  how¬ 
ever,  it  docs  not  stand  out  from  several  of  the  other  peaks. 
Similarly,  dominant  peaks  arc  not  results  obtained  for  the  cep- 
stium  of  ~4  min.  of  data  and  for  the  ccpstrum  stack  consisting 
of  9  -  25.6  second  samples. 

In  the  next  plot  down,  we  see  the  positive  effects  of  the 
stochastic  stacking  technique.  The  stochastic  ccpstrum  stack 
uses  the  identical  cepstrums  that  were  used  in  the  ccpstrum 
stack  plotted  above,  hut  by  using  this  technique  the  peak  at  the 
pP  delay  time  is  the  dominant  peak. 

In  achieving  this  the  stochastic  window  was  varied  from 
0.  to  1.6  seconds  where  the  results  plotted  arc  for  a  window 
width  of  .8  seconds.  The  results  of  the  stochastic  stacking  as 
a  function  of  the  various  time  windows  arc  plotted  in  Figure  8. 
he  see  that  the  best  results  arc  obtained  using  the  windows  in 
the  range  of  .8  and  1.2  seconds.  This  window  width  range  gave 
the  best  results  for  all  stations,  and  is  in  agreement  with  the 
maximum  delay  time  variation  expected.  That  being  the  (pP-P)  - 
(pPl  II)  diffcicncc  of  ~  1  see.  Thus,  it  may  be  possible  to 
estimate  what  window  to  use  for  a  given  event.  (Figure  A2  of 
the  appendix  shows  another  example  of  the  effect  different  sto¬ 
chastic  window  widths  have  for  the  station  FN-MO.) 

In  the  stochastic  ccpstrum  phasor  stack,  shown  in  the 
bottom  plot,  each  point  of  the  ccpstrum  was  represented  by  a 
phasor  and  the  amplitude  of  the  stochastic  ccpstrum  phasor  sum 
was  plotted.  This  technique  was  not  effective  for  this 
station  for  reasons  to  be  discussed. 


5-7 


*****•&*&* 


1 1 
I 


Window  Widths 


(seconds) 


i  i 

i  )  \ 


'"■  ■■  !\  /  „  * 

'  v  ■  V I  Wv^Va  ,. 


;  i  1 1 1 1 1 1 1 1 1  i 1  r  I.m  r  r 1 1 1 1 1 1  n  r 1  j  1 1 1 ;  1 1 ! ! , ,  „  1 1 1 1 , , , 1 1  ITT-| r i , T  m  ,  i  rj 1  rrj  1 1  nyi ii  i  r|  1 1  it|  p  1 1  (j  1 1 1 1 , 1 1 1 1,  i , ,  i , , , , , , , , , 


V  \  ./AM  f  \ 

v  v  i  vrv'-\  ,, 


/  \ 


1 1 !  1 1  m  |  i  ig 


i  >  1 1 1 1 1 1 1 1 1 1 1  ii  i  r|  1 1  u  1 1 1  r 1 1 1 1 1 1  j  1 1  it 1 1 1 1 1 1 1 1 1 1 1 1 1  ri  (i  in  1 1 1 1 1|  1 1 1 1 1 1  1 1 1 1 1  rr] i  1 1  rp  i  rrfi  mp  1 1  rp  1 1  rp  1 1 1  p  1 1 1  p  1 1  f 


VvVi\  '\  - 

v  I  V/-'V' 


\  f  N 


■  1 1  ■  i ,  1 1  ■  •  ‘  ■  1 1 1  •  1 1 1  •  •  •  •  > '  I  ■  1 1 1  •  •  1 1 1 1  ■  II 1 1 1  ■  •  r  •  •  ■  >  I  •  •  •  1 1 1 1 1  r; ! !  |  r  Tj  1 1 17 1 1  m  1 1 1  <  l  T|  1 1 1 1 1 1 1 1 1  ^  1 1 1 1|  1 1  IT|  1 1 1 1 1 1 1 1  r  1 1 1  fT[l  1 1  f  1 1 1 1 1 


/  \ 
J  \ 


. . . .  •  -  -  "  • "  ■ ' !  1 1 :  r  1 1 1 1 1 1 1  m  p  i  it|  1 1  n  p  1 1 1 1  i|i  i  r  1 1  mt  1 1 1  r1  pi  i  irp  1 1  t  1 1  h  r  1 1 1 1 1 1 1 1  ii  |  n  it  1 1 1  rr|  i  ii  r|  i  n  r;  f  riif  1 1  ir|i  •  s  r 

v\  '  I  '  I 


/  \/\ 


\ 


1 1 1 1 1 1  r  |  ii  1 1  j  i  n  ip  i  up  i  rrj  1 1 1  rp  nT|  i  rrr|  r  i  rrp  i  n  p  rnpi  it|  i  mp  nr  p  rn  1 1  |vrp  i  u  p  n  i  p  1 1  rp  1 1  rp  1 1 1  p  1 1  ip  ii  r  p  1 1  l‘p  1 1  r 


■  i  •  |>)  .•,!  |*  tJl  Cfi  ''J  IN  *-L  »■*  »-*•  h*-  i-*  i-  fO  Kj  [0  M  I 

,  •  •  -  •  •  •  •  •  •  a  r-  ro  t»i  j>  ai  a,  -'•j  to  <n  a  ■  *  fa  (• 

•  I  r,  f.  r  tv  fi  Jv  -fv  .|>  . 

i.i  m  »  rn  n,  i  ,j  4*  m  in  in  ti\  n\  tn  in  t\\  m  n»  u»  in  n\ 

rti  ifi  r .  I-  I  *  h.1  M  4.  lit  iT'i  •  I  CO  .Vi 


Stochastic  cepstrum  stack  for  various 
stochastic  window  widths  (MN-NV) 


Figure  8 
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Next  consider  the  results  for  IIN-MI]  shown  in  Figure 
Here  again  the  first  three  conventional  ccpstrums  analysis 
techniques  did  not  produce  dominant  peaks  at  either  of  the 
expected  pP  and  sP  delay  times.  In  this  ease  the  addition  of 
stochastic  stacking  w  s  not  enough  to  improve  on  the  normal 
techniques.  However,  here  the  combination  of  the  pliasor  tech¬ 
nique  with  the  stochastic  stacking  resulted  in  the  detection 
of  peaks  at  both  pP  and  sP  delay  times.  (The  phnsor  stacking 
gave  these  results  only  when  used  with  the  stochastic  window¬ 
ing.)  In  fact  the  entire  ccpstrun  pattern  of  the  stochastic 
ccpstrum  phasor  stack  is  in  agreement  with  the  ccpstrum  (bot¬ 
tom  plot)  calculated  from  a  synthetic  seismogram  for  this 
event,  Ihc  single  station  analysis,  giving  the  entire  ccpstrum 
pattern  expected  for  an  event  at  ~39  km,  is  by  itself  extremely 
strong  evidence  for  an  event  at  this  depth.  We  will  see  that 
these  results  arc  in  agreement  with  the  other  depth  determina¬ 
tions  to  be  discussed. 

The  reason  the  phasor  techniques  failed  using  data  from 
MN-NY  is  now  understood.  As  was  discussed  in  Section  2.3,  there 
is  the  possibility  of  the  phase  difference  0  to  change  by  tt 
along  the  coda.^  ^  This  was,  in  fact,  observed  in  examination 
of  the  phase  as  a  function  of  position  along  the  coda.  The 
phasor  stack  was  then  changed  to  constructively  add  phnsors  it 
out  of  phase  and  a  detection  comparable  to  the  stochastic  stack 
of  MN - N V  w a s  a c h i e v c d . 

The  amplitude  of  the  individual  ccpstrum  phasor  arrays  for 
UN-MI-  is  plotted  in  Figure  10.  Consider  the  peak  at  the  sP  de¬ 
lay  time  in  the  bottom  plot  of  the  phasor  array  sums.  By  in¬ 
specting  the  amplitudes  at  this  delay  time  for  each  of  the  9 
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Figure  10 
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component  ceps t rums ,  one  docs  not  sec  evidence  of  a  peak. 

This  is  dramatic  evidence  of  how  the  phase  information  can 
aid  in  detection;  for  in  this  case  it  is  the  consistency  of 
the  phase  of  these  phasors  and  not  their  amplitudes  which 
give  rise  to  a  dominant  peak  in  the  pliasor  summation.  This 
technique  shows  great  promise  and  should  he  investigated  fur¬ 
ther.  Modifications  such  as  increasing  the  phase  selectivity 
can  easily  be  worked  into  this  technique. 

In  figure  11  the  results  for  EN-MO  arc  plotted.  For  this 
station  there  is  no  indication  of  a  depth  phase  at  the  expected 
delay  time  from  either  the  visual  or  conventional  methods  shown 
in  Figure  11.  (In  Figure  A-6  arc  plotted  14  consecutive  25.6 
second  ccpstrums  for  UN-MO.)  Only  the  stochastic  ccpstrum 
stack  and  the  stochastic  ccpstrum  pliasor  stack  give  a  dominant 
peak  at  the  pP  delay  time. 

In  Figures  12-14  arc  results  for  stations  11V-MA,  KN-UT, 
and  JP -AT,  which  lie  in  the  Rocky  Mountains.  Depth  phase  de¬ 
tection  was  not  successful  at  these  stations  using  any  of  the 
methods . 

A  summary  of  the  results  for  the  individual  station  ana¬ 
lysis  arc  seen  in  Figures  15  and  16.  In  Figure  15  arc  the 
results  for  all  six  stations  using  methods  1  and  3.  There  is 
little  evidence  of  dominant  peaks  appearing  at  the  expected 
delay  times.  Figure  16  summarizes  the  results  of  the  stochas¬ 
tic  stacking  and  stochastic  pliasor  stacking  techniques.  Here, 
although  not  overwhelming,  is  evidence  developing  for  a  depth 
phase  detection. 
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(1st  -  25.6  Seconds) 


Cepstrum  Stack 
(9  *  25.6  Second  Samples) 


Summary  of  results  of  the  25.6  second  cepstrum 
and  the  Cepstrum  Stack 


Figure  15 
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Stochistlc  Stick 


Summary  of  Results  from  the  Stochastic  Stack  and  the 
Stochastic  Phasor  Stack 


Figure  16 


In  figure  A3  of  tlic  appendix  are  plotted  the  auto¬ 
correlations  of  the  low  pass  filtered  (cutoff  frequency  is 
equal  to  the  inverse  of  the  stochastic  time  window)  seis- 
j  ograms  for  five  of  the  stations.  To  calculate  these  auto¬ 
correlations  the  N  sampled  amplitude  values  (X^ (n) ,  (n=0 , 1 , 
N-l)  were  low  pass  filtered  using  a  running  average, 

K 

2  xM«-k> 

k=  -K 


XM0)  = 


2K+1 


where  K  is  set  at  10  such  that  the  cut  off  frequency  is 
approximately  equal  to  the  inverse  of  the  stochastic  time 
wine'  w.  N  zeroes  were  then  added  to  to  avoid  wrap  around 
problems  and  the  correlation  function  CM ( k )  for  positive 
lags  was  then  calculated  from 


PM«>  = 


ccM(k)= 


2N- 1 

1 

v 

2N 

Z 

n=0 

2N- 1 

1 

Y 

2K 

Z 

j  =  0 

NS 

I 

CCM 

M=  1 

XM  (n)  c 


2  7T  i  jn/2N 


,  j  =  0 , 1 , * . .  ,  2N 


pM(j)  e'27Tijk/2N,  k=0 , 1 , .  .  .  ,N- 1 


M(k)  -  «.„(*)  '  WM  .  k=0>1 . 
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In  Figure  A4  of  the  .appendix  arc  plotted  the  cep- 
strum  stack  of  the  low  pass  filtered  seismograms,  where 
a  pP  detection  similar  to  the  one  obtained  by  the  stochastic 
techniques  is  apparent  for  MN-NV.  This  is  the  only  ease  in 
which  one  of  the  conventional  techniques  gave  rise  to  a 
dominant  peak  at  an  expected  delay  time.  For  this  case 
low  pass  filtering  had  much  the  same  effect  as  the  stochastic 
windowing.  Figure  A5  of  the  appendix  shows  the  combined 
best  results  of  methods  4  and  5.  Also  tried  but  not  plotted 
was  analysis  using  low  pass  filtered  data  with  the  phasor 
technique  without  using  stochastic  stacking.  This  method 
was  not  successful  for  any  of  the  six  stations. 


Summarizing  the  results  of  all  methods  discussed  so  far, 
depth  phase  detection  was  achieved  at  three  of  the  six  stations 
(one  of  which  gave  both  pP  and  sP)  using  the  new  techniques, 
whereas  the  conventional  methods  achieved  pP  detection  at  one 

station . 


3.2  Pc pth  Phase  Pctcction  Using  Combined  Data  from  Sevcial 
Stations 

Next  we  show  how  further  improvements  in  depth  phase  de¬ 
tection  can  be  achieved  by  using  the  combined  data  of  several 

stations. 


Return  to  Figure  16  where  one  observes  that  for  sta¬ 
tions  1IV-MA ,  MN-NV,  KN-NT,  pP  is  detected  at  only  MN-NV  and 
sP  is  not  detected  at  any  of  these  stations  using  stochas- 

■*  i  .  a  _  1 


t  i  c  s 


•;  t  nek  i  lUl 


For  each  of  these  stations  we  calculated 
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an  equal  number  of  ce  strums  from  consecutive  25. 0  second 
portions  of  the  seismograms.  In  figure  17  arc  the  results 
of  stochastically  averaging  these  ccpstrunis  of  these  sta¬ 
tions  using  different  coda  lengths  (the  trend  has  been  re¬ 
moved  from  the  averaged  ccpstrum  plotted).  These  three  sta¬ 
tions  arc  at  approximately  the  same  distance  from  the  event 
(^45°).  We  note  a  clear  detection  of  b(  th  the  pP  and  sP  delay 
times  when  5.0  minutes  of  the  codas  is  used.  Also,  note  that 
the  detection  of  pP  and  sP  improves  as  more  of  the  information 
contained  in  the  coda  is  used.  This  is  evidence  that  these  sto¬ 
chastic  techniques  arc  making  constructive  use  of  the  additional 
seismic  dcptli  phase  information  contained  in  the  coda.  (Similar 
results  were  obtained  by  stochastically  stacking  the  ccpstrum 
phasors  of  the  three  stations.) 

The  results  can  be  further  improved  with  the  addition  of 
the  seismic  data  recorded  at  the  stations  EN-MO,  and  IiN-ME.  In 
Figure  19  arc  th c  results  of  averaging  data  from  the  five  sta¬ 
tions  I1V-MA ,  MN'-N’V,  KN-IJT,  EN-MO  and  UN’ -ME  using  coda  lengths 
-2.5  and  -5.0  min.  W’c  obtain  a  clear  detection  of  pP  using 
2.5  min.  of  the  coda  and  by  using  5.0  min.  the  sP  detection  be¬ 
comes  clear.  But  note  that  the  two  additional  stations  arc  at 
62°  and  66°  compared  with  -45°  for  the  other  three.  The  ccpstral 
peak  position  would  th ere  fore  be  expected  to  differ  by  ~ . 5  sec¬ 
onds  between  the  two  sets  of  stations.  One  can  adjust  for  these 
station  distance  differences  by  shifting  the  ccpstrums  for  EN-MO 
and  I  IN -ME  by  --.5  seconds  before  stacking.  A  further  improve¬ 
ment  in  clcarity  in  the  detection  of  pP  and  sP  delay  times  is 
then  observed  from  Figure  19. 
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Stochastic  cepstrum  stack  using  25.6  seconds  samples 
of  the  3  stations  HV-MA,  MN-NV,  KN-UT 


Figure  17 


Fhc  results  shown  in  Figure  18  are  very  encouraging 
since  we  have  obtained  an  excellent  detection  of  the  pP 
and  a  good  detection  of  the  sP  delay  times  for  an  event  in 
which  both  the  visual  and  the  standard  techniques,  which 
we  investigated,  gave  vc.y  little  indication  of  the  presence 
of  these  phases.  These  results  also  give  promise  of  signi¬ 
ficant  improvements  in  depth  phase  detection  for  small  events 
h a \ i n g  poor  signal  to  noise  ratios  when  these  events  arc  re¬ 
corded  at  several  stations.  The  stochastic  averaging  allows 
one  to  make  constructive  use  of  information  contained  in  the 
codas  recorded  at  different  stations  as  well  as  different 
portions  of  a  given  coda. 
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Stochastic  cepstrum  stack  using  25.6  seconds  samples 
of  the  5  stations  HV-MA,  MN-I1V,  KN-IJT ,  CN-MO,  HN-MH 


Figure  18 
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CONCLUSION 


This  research  effort  has  been  directed  at  determining 
whether  certain  new  statistical  signal  processing  techniques, 
applied  to  the  entire  seismogram,  could  significantly  enhance 
seismic  depth  phase  detection.  The  new  analysis  techniques 
allow  one  to  make  constructive  use  of  the  additional  depth 
Phase  information  contained  in  the  coda,  rather  than  restrict¬ 
ing  the  analysis  to  the  first  arrival  portion  of  the  seismogram. 
Results  of  the  analysis  of  a  single  event  at  six  stations  indi¬ 
cates  that  major  improvements  in  depth  phase  detection  arc  possi 
blc . 


Two  new  techniques  were  introduced  in  the  depth  phase  ana¬ 
lysis.  The  first  of  these  techniques  allows  for  depth  phase 
del  a)’  time  variations  associated  with  the  later  seismic  arrivals 
comprising  the  coda.  The  second  of  these  techniques  utilizes 
the  phase  consistency  between  the  direct  and  surface  reflected 
sc  i sm i c  arri val s . 

Lor  the  event  analyzed  the  depth  phases  were  not  visually 
identifiable  from  the  seismogram.  Results  of  the  ccpstral  ana¬ 
lysis  incorporating  these  two  new  techniques  gave  rise  to  domi¬ 
nant  ccpstral  peaks,  for  three  of  the  six  stations,  which  were 
consistent  with  a  39  km  deep  event.  At  one  of  the  stations  the 
entire  ccpstral  pattern  expected  for  a  seismogram  containing  the 
pi1  and  sP  arrivals  was  obtained.  These  results  w'erc  compared 
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i t h  the  analyses  using  several  conventional  depth  phase 
detection  techniques  including  the  cepstrum,  stacked  cep- 
strums,  auto- correlations  and  low  pass  filtering  the 
seismograms  before  performing  these  analyses.  One  of 
the  conventional  techniques  gave  rise  to  an  indication 
of  the  pP  delay  time  at  one  of  the  stations  similar  to 
that  obtained  by  the  new  techniques. 

Results  of  stochastically  averaging  the  data  from 
several  stations  demonstrated  detection  of  depth  phases 
which  were  not  apparent  from  analysis  at  any  of  the  in¬ 
dividual  stations.  These  results  also  clearly  indicated 
how  the  dcptli  phase  detection  improved  as  more  of  the 
coda  was  used  in  the  analysis. 

The  new  dcptli  phase  detection  techniques  have  been 
thus  far  applied  to  only  a  single  event.  The  analysis 
gave  clear  indications  of  pP  and  sP  delay  times  consis¬ 
tent  witli  a  39  km  deep  event  and  represented  a  major 
improvement  over  conventional  techniques  for  this  case, 
however,  the  depth  of  this  event  was  not  considered  to 
have  been  conclusively  established.  The  next  research 
effort  should  therefore  lie  to  apply  these  new  techniques 
to  a  suite  of  events  having  well  establish'd  depths  and  to 
further  develop  these  techniques  using  this  data. 
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Demonstration  nf  stochastic  stacking 
using  low  pas  filtered  synthetic  data 


S-o{ 


Figure  A-l 
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Stochastic  cepstrum  phasor  stack  for  various 
stochastic  window  widths  (EN-MO) 


Figure  A-2 
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Auto-Correlation  Stack  of  Low  Pass  Filtered  Data 


(9-25.6  second  samples) 


Figure  A-3 
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Cepstrum  Stack  of  Low  Pass  Filtered  Data 
(9  -  25.6  second  samples) 


Figure  A-4 
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Combined  best  results  of  the  stochastic  cepstrum  stack 
and  stochastic  cepstrum  phasor  stack  techniques 
(9  -  23.6  second  samples) 


Figure  A-S 
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